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ABSTRACT 

The discovery of SN 2002ic and subsequent spectroscopic studies have led to 
the surprising finding that SN 2002ic is a Type la supernova with strong ejecta- 
circumstellar interaction. Here we show that nearly 1 year after the explosion 
the supernova has become fainter overall, but the Ha emission has brightened 
and broadened dramatically compared to earlier observations. We have obtained 
spectropolarimetry data which show that the hydrogen-rich matter is highly as- 
pherically distributed. These observations suggest that the supernova exploded 
inside a dense, clumpy, disk-like circumstellar environment. 

Subject headings: stars: individual (SN 2002ic) - stars: supernovae - stars: spec- 
troscopy - stars: polarimetry 



1. Introduction 

SN 2002ic (Hamuy et al. 2003; Wood-Vasey et al. 2002; Wood-Vasey et al. 2003) is an im- 
portant supernova (SN). Spectroscopy shows that it is a Type la SN with strong hydrogen 
emission arising from the ejecta-circumstellar interaction. Hamuy et al. (2003) point out 
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that SN 2002ic bears a strong spectroscopic similarity to SN 1997cy, previously identified as 
a Type Iln (Germany et al. 2000; Turatto et al. 2000). 

We present in this Letter spectropolarimetry observations of SN 2002ic obtained nearly a 
year after the explosion. This is the latest epoch for which spectropolarimetry has been 
obtained for any SN. We discuss the nature of the polarization and the structure of the 
circumstellar environment surrounding SN 2002ic. 

2. Observations 

We observed SN 2002ic on 2003 July 8.4 (UT), July 29.4 (UT), Jul 31.3 (UT), and Sept 29 
(UT) using the Very Large Telescope (VLT) of the European Southern Observatory (ESO) 
with the FORSl instrument in spectropolarimetry mode. The conditions during the July 

8.4 observations were photometric, and the seeing was around 0.7." The slit width used in 
these observation was 1." Four separate exposures each of 30 minutes were obtained with 
the retardcr wavcplate at position angles of 0, 45, 22.5, and 67.5 degrees. The July 8.4 
(UT) data provide both reliable spectropolarimetry and spectrophotometry. The conditions 
for July 29.4 (UT) and 31.3 (UT) were less ideal. Only exposures at waveplate position 
angles of 0, 45, and 22.5 degrees were obtained on July 29.4 (UT), but all four position 
angles were again observed on July 31.3 (UT). The July 29.4 (UT) and July 31.3 (UT) 
data were combined to produce the spectropolarimetric results. The flux calibrations of the 
July 29.4 (UT), 31.3 (UT) and Sept. 28 (UT) data are less reliable but not by more than 
10% as found from the scatter of the final reduced data set of each individual exposure. 
The spectral profiles show little evolution during these observations. Figure 1 shows the 
spectroscopic data derived from the July 8.4 (UT) observations together with the spectra 
of SN 1997cy (Turatto et al. 2000) and SN 1999E (Rigon el al. 2003) at comparable dates. 
The most prominent spectral line is the broad sharply-peaked Ho; emission that is typical 
of a Type Iln SN. SN 2002ic remains strikingly similar to SN 1997cy and SN 1999E at later 
epochs (Hamuy et al. 2003). 

2.1. Time Evolution and Comparisons to Other Type Iln Supernovae 

The HcK line in our observations is broader than measured in earlier observations (Hamuy et al. 2003). 
The full width at zero intensity (FWZI) can be traced to about 20,000 km s~^ compared 
to an earlier value of about 4000 km s~^. The Ha line profile shows a strong central peak 
and is significantly different from that of the broad Mg I] 457.1 nm line. Several hues from 
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neutral helium and doubly ionized oxygen were also detected. The FWZI of these lines are 
typically around 5,000 km s~^. The profiles of these lines are actually very similar to the 
central peak of the Ha line. The broad wings, if present, may not be detectable because of 
line blending and the faintness of these lines. 

The flux of the Ha line increased by a factor of more than 6 compared to that of 2002 
Nov 11 to 2003 Jan 9. The light curve of the Ha line is shown in Figure 2, together with 
those of SN 1997cy (Turatto et al. 2000) and SN 1999E (Rigon el al. 2003) for comparison. 
The Ha fiuxes derived from our spectra are 2.2, 1.9, 1.8 and 1.7 in units of xlO^^'' ergs 
cm-2 s-i for observations of July 8.4 (UT), July 29.4 (UT), July 31.3 (UT), and Sept 28 
(UT), respectively. The brightening of the Ha line after optical maximum is not unique to 
SN 2002ic, but is also observed for similar events such as SN 1988Z (Turatto et al. 1993), 
SN 1997cy (Turatto et al. 2000), and SN 1999E (Rigon el al. 2003). The spectroscopic data 
show no remarkable evolution during these observations. 

2.2. Spectropoleirimetry 

To investigate the details of the polarization spectral variation, it is necessary to identify 
the dominant axis of asymmetry in the plot of Stokes parameters, Q versus U, and to 
decompose the data into components along the dominant axis (P^, hereafter) and the axis 
orthogonal to the dominant axis, or the secondary axis (Pq, hereafter) (Wang et al. 2001; 
Wang et al. 2003). The advantage of identifying the dominant axis is that the spectral 
features of P^ will not be sensitive to the uncertainties of interstellar polarization (ISP) if 
the majority of the polarization is produced by a global asphericity of the emitting region. A 
weighted fit of the data points on the Q-U plot binned to 10 nm shows that the dominant 
axis is at a position angle of —10° ± 2. 

The polarization due to Galactic dust particles in the field of SN 2002ic should be less than 
0.15% (Heiles 2000) and cannot be the cause of the polarization. Dust scattering in the 
host galaxy may contribute. To determine the amount of host interstellar dust polarization, 
we assume that the observed polarization of the central, narrow peak of the Ha line (from 
6553 - 6573 A) is due completely to interstellar dust. This assumption is justified because 
line photons seen in the H-peak are produced by recombination and, thus, are intrinsically 
unpolarized (in the absence of large magnetic fields). In addition, the narrow-line photons, 
unlike the broad-lines photons, are not polarized by Thosmson scattering in the circumstellar 
disk because they emerge predominantly from regions with negligible electron-scattering (the 
thermal velocity of electrons is sufficiently high that the scattered photons are shifted to 
the broad wings of the Ha profile). We thus take the ISP to be that corresponding to 
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the central 20 Aof the Ho; emission peak. This gives an estimate of the ISP of about 1.6 
percent with components Q = 0.9% and U = —1.3%. With this ISP we determine the net 
polarization to be about 1 percent for wavelength regions outside the narrow Ha emission 
peak. The dominant and secondary polarization components after correcting for ISP are 
shown in Figure 3. It can be seen that across the Ha line center the dominant component 
Pd is indeed different than that of other wavelength regions. 

The polarization is different from what was found in other SN la such as SN2001el (Wang et al. 2003) 
which shows a strong feature at 615 nm corresponding to the strong line of Si II. The polar- 
ization in SN 2002ic is likely to be related to the circumstellar medium (CSM) rather than 
to the explosion process itself. 

3. The Circumstellar Matter Disk 

3.1. Constraints from Ha Emission Line 

At the dates of these observations, SN 2002ic is powered by the high energy photons from the 
ejecta-CSM shock. The velocity of the ejecta-CSM shock is found to be around 11,500 km 
s~^ from the width of the emission lines associated with the ejecta, such as the Mg I] 457.1 
nm line and the Ca II IR triplet. As noted above, the Ho; hne has a remarkably different 
profile than the Mg I] line, suggesting different origins of these emission lines. Also, the He 
I and the [O HI] lines are typically narrower than the Mgl] line with profiles similar to the 
central peak of the Ha. This is important evidence that the narrow component of the Ha 
and the He I lines are produced in the CSM and that the composition of the SN ejecta, as 
reflected in the Mg and Ca lines, is consistent with that of a typical Type la. 

If the CSM is due to a stellar wind, the luminosity of the Ha from the unshocked wind is 
~ 4.4 X lO'^^rsh ig^^M^gUj"^ erg s~^, where Tsh,w is the radius of the shock in units of lO"*^^ cm 
and we have taken the emissivity for Case B recombination at a temperature of ~ 10,000 
K. The emissivity is a weak function of temperature and changes slightly with Ha opacity 
(Xu et al. 1992). This is a decreasing function of time and cannot explain the rise of the Ha 
flux. If the CSM were dense enough, the shocked CSM may become radiative and produce 
enough Ha emission to account for the observed flux; however, the line proflle from a thin 
shell of material moving at the shock speed of 11,500 km s~^ would be flat-topped and hence 
inconsistent with the observed Ha line profile. The mass of hydrogen required to produce 
the observed Ha luminosity is ~ 6M0 YoSchF^ • ^'^'^^ realistic modeling of Ha emission may 
alter this estimate but is unlikely to change the qualitative conclusion that a few solar masses 
of dense hydrogen-rich matter is required to explain the Ha luminosity. 
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3.2. Constraints of the Spectropolarimetry 

The configuration we encounter here is different from those for normal SN la where the 
polarization is produced by electron scattering from an aspherical photosphere (Hoflich 1991). 
Mass loss from AGB stars is known to be highly asymmetric and predominantly on the 
equatorial plane. The polarization is due to scattering by the dense CSM of photons originate 
from the SN at the center. Such systems have been studied in great detail in models of 
polarizations of Be stars. In these models the opening angle of the scattering CSM has the 
shape of disk with openning angles typically of the order a few degrees in order to repoduce 
an observed polarization of order 0.5 - 1% (McDavid 2001; Melgarejo el al. 2001). Such 
a geometry is in good agreement with the required several solar masses of hydrogen with 
densities above 10^ cm~^ that contributes to the Ha emission. The dense CSM must have a 
small volume filling factor. 

4. Ejecta-CSM Interaction and the Ha Line Profile 

If the mass loss during the AGB phase were predominantly in the equatorial plane, we 
would expect the shock of the ejecta-CSM interaction to be hour-glass shaped. The disk of 
unshocked CSM could be ionized by the hard photons from the ejecta-CSM interaction and 
give rise to the central peak of the Ha emission whereas the shocked CSM could be radiative 
and give rise to the broad Ha wings. Detailed calculations are needed to see if this scenario 
can account for the observed line profile and the increase of Ha luminosity while the overall 
luminosity decreases. 

Another possibility is that broad wings of Ha are formed by scattering, rather than emission 
from shocked matter. If the CSM is optically thick to electron scattering, electron scattering 
can effectively broaden the Ha line and produce a center-peaked profile (Chugai 2003). An 
example of such a profile is shown in Figure 2. For this simple model spherical symmetry 
and a temperature of 30,000 K, presumably produced by the fiash ionization of the CSM, 
are assumed. Note that this model can produce the smooth line profile with both the narrow 
peak and the broad wings whereas such an agreement would require careful tuning of a model 
where the broad wings were shock-accelerated material and the narrow core were from nearly 
stationary CSM. In this model, photons that remain within the core have not been scattered 
and is thus unpolarized. 

It is likely that the flattened CSM is not uniform, but clumpy so that the Ha line is dominated 
by emission from the densest clumps. The clumps may have densities of about 10^ — 10^ 
cm~^ with sizes of about 5 x 10^^ cm at which point they would be marginally optically thick 
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to Thompson scattering and hence provide the sort of hne profiles illustrated by the simple 
model in Figure 2. As noted above, about 1 Mq of matter is needed in the clumps/disk to 
reproduce the observed Ho; luminosity. 

If the disk and dense clumps extend to ~ 200 hght days, ~ 3 x 10^^ cm, the rising flux of Ha 
could be the result of the time delay effect. At early times, only a small fraction of the dense 
clumps would be visible and only at around day 200 when these observations were made 
would the majority of the CSM clumps be observed, therefore giving rise to the increased 
Ha luminosity. In the mean time, the emission from the ejecta-CSM interaction should show 
a constant decrease of luminosity after reaching maximum a few weeks past explosion, as 
observed. 



5. Conclusions 

The strong growth in the Ha luminosity, the Ha line profile, and the polarization we observe 
for SN 2002ic are consistent with several solar masses of clumpy material arrayed to distances 
of ~ 3 X 10^'' cm in a manner that is globally asymmetric. The CSM may be related 
to the dense disks or rings observed in some proto-planetary nebulae (Zijlstra et al. 2001; 
2; Su et al. 1998). The mass loss was found to be highly asymmetric at this stage. For 
example, IRAS 16342-3814 (Zijlstra et al. 2001) shows a dense torus with mass 0.1 Mq and 
density of 10^ cm~^. Alternatively, a similar configuration may be produced during the 
binary evolution leading to the thermonuclear explosion. During a phase of rapid expansion 
of the donor star very high mass loss may be channeled through the outer Langrange points. 
Wind interaction (Kwok 1993) can create dense shells or rings that may naturally explain 
the observed structure (Kwok 1993; Young et al. 2003). 

The star that exploded as SN 2002ic could be either the dying post AGB star (Swartz et al. 1991) 
or a pre-existing white dwarf companion around it. The large radial extent of the disk seems 
to be incompatible with a model in which the hydrogen represents an expelled common en- 
velope at a distance of ~ 10^^ cm (Livio & Riess 2003). The proto-planetary nebula phase 
lasts only for a few hundred years. This may account for the rarity of events like SN 2002ic, 
but raises questions of how the SN is timed to precisely this phase. Alternatively, a struc- 
ture with similar properties may be produced during the binary evolution leading to the 
explosion, although questions of the phasing remain. In any case, events such as SN 2002ic 
suggest that SN la can occur in a wide range of CSM environments. 

We are grateful to ESO for the generous allocation of observing times. We are indebted to 
Massimo Turatto for providing data of SN 1997cf and SN 1999E in digital format. 
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Fig. 1. — The spectra of SN 2002ic are compared with those of SN 1997cy and SN 1999E. 
From top to bottom, the spectra are: SN 1999E, SN 2002ic on 2003 July 8, Jul. 31, and 
Sept 28 (UT), and SN 1997cy. The SN 1997cy and SN 1999E spectra are at epochs where 
the Ha line is close to maximum (see Figure 2). 
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Fig. 2. — (a) The Ha luminosity of SN 2002ic compared with those of SN 1997cy and SN 
1999E. The sohd hne illustrates a linear increase with time of the line flux. SN 2002ic shows 
comparable Ha emission to SN 1997cy and SN 1999E. The dates of explosion for SN 1997cy 
and SN 1999E are uncertain and the data points can shift by more than 20 days on the 
horizontal axis, (b) The Ha line profile (solid line) of SN 2002ic is compared with the H/3 
profile (dotted line). Also shown is a line profile calculated for a spherical H II region with 
electron temperature of 30,000 K, and Thomson scattering optical depth of 4 (dashed line). 
The overall agreement of the line profile suggests electron scattering plays a major role in 
forming the line profiles. 
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Fig. 3. — The polarization vectors (middle and lower panel) of SN 2002ic are shown together 
with the flux spectrum (top panel). The Stokes parameters are projected onto the dominant 
axis (Prf) and the axis orthogonal to the dominant axis (Po)- The thick solid line shows the 
combined data of all observations, the thin red and green lines show the July 8.4 (UT) and 
July 31.4 (UT) observations, respectively. The dominant axis is constructed by a weighted 
linear fit of the data points on the Q-U plot, and the position angle of the dominant axis was 
found to be 10o±0.2. In constructing and Pq, we have also assumed that the central peak 
of the Ha line is unpolarized to deduce the ISP component. The ISP used in calculating Pa 
and Po are Q = 0.9% and U = —1.3%. The long dashed lines illustrate the mean values 
of the polarization vectors, and the short dashed line shows the location of the Ha line. 
The Ha line clearly shows a different degree of polarization when compared to the overall 
polarization of the wavelength region from 300 nm to 600 nm. 
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